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Abstract 

Given an quantum dynamical semigroup expressed as an exponential superoperator acting on a 
space of ^-dimensional density operators, eigenvalue methods are presented by which canonical 
Kraus and Lindblad operator sum representations can be computed. These methods provide a 
mathematical basis on which to develop novel algorithms for quantum process tomography, the 
statistical estimation of superoperators and their generators, from a wide variety of experimental 
data. Theoretical arguments and numerical simulations are presented which imply that these 
algorithms will be quite robust in the presence of random errors in the data. 
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I. INTRODUCTION 



The statistical estimation of superoperators from experimental data is variously known as 
"quantum channel identification" |l] or "quantum process tomography" (QPT) 0. While 
this task is important throughout experimental quantum physics, it is an essential component 
of on-going efforts to develop devices capable of reliable quantum information processing and 
transmission. At the same time, it is only through these efforts that it is now becoming pos- 
sible to observe and control quantum systems with the precision needed to collect sufficient 
data for QPT. At the time of writing, however, very few experimental efforts to systemat- 
ically determine the complete superoperators of natural or engineered quantum processes 
have been carried out. An instructive example may be found in ||, where the QPT proce- 
dure detailed in || was applied to NMR data on the two-qubit molecule chloroform. This 
was followed by fitting a specific decoherence model to the superoperators thereby obtained 
at multiple time points, in order to estimate the decoherence rates in the model. 

The goal of the present paper is to give a reasonably complete and self-contained account 
of the mathematics needed for robust QPT, assuming for the most part that the quan- 
tum dynamics may be aptly modelled as a quantum dynamical semigroup (QDS). A QDS 
describes the evolution of a general open quantum system under the Born-Markov approxi- 
mations 0, |5|, and as such is sufficient to cover most of the systems currently being used 
or developed for quantum information processing and transmission. By "robust" , we mean 
that the QPT results will not be sensitive to random errors in the data, which is critical 
since these data are often difficult to obtain and significantly contaminated by noise and 
other errors. In addition, it is desirable to avoid model fitting and instead to determine the 
complete superoperator making no prior assumptions about it, although this significantly 
increases the number of parameters to be estimated. 

The robustness of our approach is obtained primarily by using the orthogonal projection 
of an arbitrary Hermiticity-preserving superoperator or QDS generator onto the convex cone 
of completely positive superoperators and their generators ||, [|, @- Of necessity, therefore, 
this account will rederive much that is already known about quantum dynamical semigroups 
as well as more general completely positive superoperators, using a consistent notation, fixed 
operator basis, and a standard set of matrix tools []7|, §]. These derivations do not involve 
qualitative physical arguments (coarse-graining, separation of time scales, etc.), but only 
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the mathematical definitions of the objects involved, and extend much of our earlier work 
on the "Hadamard" representation, which exists for any "diagonal" superoperator, to more 
general completely positive superoperators and QDS generators ||. 

The main results will be eigenvalue methods by which the projection of an arbitrary 
Hermiticity-preserving superoperator or QDS generator onto the convex cone of completely 
positive superoperators can be computed. These projections will be shown to yield certain 
canonical Kraus and Lindblad representations of completely positive superoperators and 
QDS generators, respectively, which may be novel and are certainly not well-known. The 
explicit form of the involution which identifies a Hermiticity-preserving superoperator with 
a quadratic form (or Hermitian supermatrix) , herein denoted by "C/ioi", also appears to 
be new (see Corollary |2] ff.). It is not the intention of this paper to give a single fixed recipe 
for QPT, because any such recipe must depend to some extent on the nature of the data to 
be analyzed. Nevertheless, a simple example will be given using simulated data plus added 
random noise, which should make it clear how such recipes can be derived from these results 
and further demonstrates that such recipes may be expected to be robust. 



II. BACKGROUND ON QUANTUM DYNAMICAL SEMIGROUPS AND THEIR 
REPRESENTATIONS 



This section provides the essential background on quantum dynamical semigroups needed 
in the remainder of the paper, and in addition defines the basic mathematical operations and 
notation to be used throughout the paper. A quantum dynamical semigroup (QDS) ||, 0, |6| 
constitutes a bounded one-parameter family of "superoperators" S = S(t; • ) acting linearly 
on a space of self-adjoint "density" operators p, and satisfying S( t+t' ; p ) = S( t ; p)S( t' ; p) 
for all p and t, t' > 0. Assuming that p acts in turn on a complex Hilbert space of dimension 
iV < oo, a general means of representing a QDS is as a Kraus operator sum ||10| , namely 

M 

pit) = S(t;p) = J2 S m(t)pSUt) , (1) 

m=0 

where one may take M < N 2 , the S m act on the same Hilbert space as p = p(0), and 
the dagger (f) denotes the adjoint. This ensures not merely that S preserves the positive 
semidefiniteness of the density operator p, but moreover that it is completely positive, mean- 
ing that the trace over any other quantum system on which S acts trivially is again a positive 



semidefinite operator, as expected for any physically realizable process [loc. cit.]. 

On identifying the Kraus operators S m with a matrix representation thereof S m , a well- 
known result regarding Kronecker matrix (or tensor) products M, Q implies 

M 

col{p{t)) = S(t)col(p) = {S m (t)®S m (t))col(p) , (2) 

m=0 

where col(p) denotes the result of stacking the columns of the corresponding density matrix 
p in left-to-right order on top of one another to get a single column vector of dimension N 2 , 
the overline denotes the complex conjugate, "<S>" the Kronecker product and juxtaposition 
denotes matrix multiplication. Although this result (which can be proved by straightforward 
index gymnastics) is often neglected in theoretical treatises on open quantum systems, it 
is extremely useful for computational purposes. In particular, it converts the two-sided 
operations in the Kraus operator sum to one-sided matrix operations, thereby providing 
a matrix representation of the one-parameter semigroup S. It further makes clear that a 
completely general linear transformation T of the "Liouville" (matrix) space <C NxN can also 
be written in operator sum form as 

N 2 -l 

X' = T(X) = r ™ T ™ X T « > ( 3 ) 

m,n=0 

where X G C NxN , the T m are a matrix basis thereof, and r mn G C are coefficients, since 

/ n 2 -i \ 
coJ(X') = I r mn T n ®T m I col(X) = TcoZ(X), (4) 

\ m,n=0 / 

and [T n ® T m \ < m, n < N 2 — 1] constitutes an induced basis for the space of "super- 
matrices" C N * N . Clearly, T preserves Hermiticity if and only if the matrix of coefficients 
b~mn]m,n=o is Hermitian. 

The semigroup property S(t + t'; p) = S(t; p)S(t'; p) is of course not assured by the 
existence of a Kraus operator sum representation, but it is equivalent to the existence of 
a constant superoperator Q G C n2xn2 such that S(t) = Exp(-Qt) for all t > 0, where 
ll Exp(—<3 1) = 1 — Q t + \Q 2 1 2 + ■ • ■ " is the corresponding exponential superoperator (see 
Ref. fllj and citations therein). In general, however, such an exponential will not possess 
a Kraus operator sum representation, even if the real parts of the eigenvalues of Q are 
nonnegative (ensuring that the evolution is bounded). General necessary and sufficient 
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conditions for a bounded one-parameter family of superoperators to be a QDS were first 



derived independently by Lindblad [|12|] and by Gorini, Kossakowski and Sudarshan |T3 
and require that the derivative p can be written in the so-called Lindblad form, 

M 

p(t) = C(p) = i[p(t),H] + ( L mP(t) L l ~ \LlL m p(t) - \p(t)LlL m ) , (5) 

m=0 

where M < N 2 as above, and both the Hamiltonian H and Lindblad operators L m are 
time-independent. The superoperator C itself is called the Lindbladian. Translated into 
matrices, this implies that the decoherent part Q of the Lindbladian can be written as 

M 

(L m ®L m - |I®(LLL m ) - i(L^L m )®l) , (6) 



m=0 



where I is the N x N identity matrix. 



III. A CANONICAL KRAUS OPERATOR SUM REPRESENTATION 

Although superoperators on Liouville space can be represented with respect to an arbi- 
trary supermatrix basis, as in Eq. (|3]), any Liouville space basis induced by an arbitrary 
Hilbert space basis can be regarded as the basis of elementary matrices E.^ (with a "1" in 
the ij-th position and zeros elsewhere), which has the advantage identifying the coefficients 
and the supermatrix elements. For example, one can write the transpose of an arbitrary 
N x N matrix in operator sum form as 

N-l N-l 
i,j=0 i,j=0 

or equivalently, as 

N-l 

col(X T ) = ( ^ Eji® Ey)cof(X) ee Kcol(X) . (8) 

i,j=0 

The supermatrix JC plays a important role in what follows. It is easily seen to be both 
symmetric and orthogonal, i.e. involutory. Using the relation E^- = e^ej (where e i; are 
the elementary unit column vectors) together with the mixed product formula (A®B)(Ccg> 
D) = AC ® BD U|, we can also show that it has the interesting property of swapping the 
order of the factors in a Kronecker product, since 

N-l N-l 

7C(X®Y)7C = EjiSEy) (X®Y) E ft ® E w ) 

i,j=0 ' k/=0 
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AT-1 

(E i< XE tt )®(E y YE w ) 

2V-1 

E (e,(e7X e,) ej) ® (e^ejY e k ) ej) 

i,j,k,e=o 

N-l N-l 



(9) 



( Y ik E ifc ) ® ( X « E « 

j,fc=0 i,£=0 



Y®X 



where the matrix elements have been denoted by = e^Xe^ and Yjk = ejYefc. 

We now use the relation coZ(xy T ) = y <g> x for arbitrary column vectors x, y to show 
how the matrix JC also gives us the col of a Kronecker product of matrices as a Kronecker 
product of their respective co/'s. 

Lemma 1 Given any two N x N matrices X, Y, we have 

col(X) (g) col(Y) = (l®/C®l)coZ(X®Y) , (10) 
where K is defined as in Eq. (|5[). 

Proof. Applying the definitions, we obtain: 

N-l 



(l<8> 7C® I) coZ(X® Y) = ((I® E^) ® (E y ®l)) coZ(X® Y) 

/AT-1 N 
= col I ^ (Ey<g>l) (X®Y) (l<g>Ey) 



AT-1 



= C0llj2 ((EyX)® (YEy)) (11) 

\ i,i=0 / 

(AT-1 

^ ((e,(eJX))®((Ye 4 )e J 1 

= co ' (( E ( e > ® ( Y e *))) ( E x ) ® e J)) 

= co/(co/(Y) co/ T (X)) = coZ(X) ® coZ(Y) 

QED 

Corollary 2 Wii/i everything defined as in the Lemma, 

N-l 

col(Y)col\X) = ^ (E ii ®l)(X®Y)(l®E ii ) . (12) 

i,j=0 
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Proof. Just apply the inverse of the col operation to the second and last lines of the proof 
of the Lemma, and add a complex conjugation to account for our use of "t" instead of "T". 
QED 

The "super-superoperator" on the right-hand side of Eq. flT2p maps any iV 2 x iV 2 su- 
permatrix <S, acting on iV x jV matrices X as Scol(K), to a new supermatrix T = 
X]fj=o(Eij <S> I) S (I ® Ey), the elements of which are a permutation of those of S. The 
Corollary shows that if S is a sum of Kronecker products, as in Eq. (|]), then T is a sum of 
the corresponding rank one dyadic products, as in Eq. (jl2|). Thus, while Eq. (Q) allows us 
to construct a supermatrix representation from an operator sum, we are now able to give a 
procedure for going in the other direction. 



Proposition 3 Let S,T e C N xN with T = J2i,j=o ® I ) S i 1 ® E *i)> and let 

N 2 -l 

r = vnW = > ( 13 ) 

fc=0 

be the singular value decomposition of T (where v^,w^ are the columns of the unitary 
supermatrices V, W, respectively, and LOf. > are the singular values). Then for any X G 

(N 2 -l \ 
X WfcVfcXWN , (14) 
k=0 ) 

where coZ(V^) = ffc, coZ(Wfc) = and the symbol should be read as "applied to". 

Proof. This follows immediately from our foregoing observations, together with the fact 
that the super-superoperator is involutory, since 

N-l / N-l \ 

X (E w ® I) I X) ® I) ^ (I ® Ey) (I <g> E M ) 

JV-1 

= X Su8jk(Ekj®l)S(l®Eu) (15) 

i,j,k,e=o 

,N-1 v , JV-1 v 

\ j=0 ' ^ i=0 



QED 
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The matrices {V^} and are not generally unitary, but each of these two sets forms a 

basis for C NxN , and each is orthonormal with respect to the Hilbert-Schmidt (or Frobenius) 
inner product (X, Y) = tr(X^Y). By expanding the right-matrices as linear combina- 
tions of the left Vfc, one can rewrite the action of T on X in the more symmetric form given 
in Eq. ([|). Thus we have obtained a general means of converting a supermatrix S acting on 
columnized matrices coZ(X) to operator sum form. Much of the foregoing can of course be 
extended to nonsquare matrices C MxN , but we shall have no need of that here. 

In the case that <S is a (super)matrix representation of a quantum dynamical semigroup, 
the matrix T derived from it has considerably more structure, as we shall now show. 

Proposition 4 With everything defined as in Proposition [| ; the derived supermatrix T can 
he written as 



N-l 

col{Sij) cof (Ei 

i,j=0 



N-l 
i,j=0 



S(En) <S(E 12 ) 
<S(E 2 i) 5(E 2 2; 



(16) 



where «S(Ey) G C NxN is defined by coZ(«S(Ejj)) = ScoZ(Ejj), and is the ij-th N x N 
block of the supermatrix S . 

Proof. The first equality in Eq. (|T6| ) follows immediately from Corollary |2| together with 
the obvious fact that «S = i=o ® ■ ^° P rove the second, we first note that for any 
< k,£ < N, 



Scol(E 



k() 



N-l 



Y Eij ® Sij )col(E M ) 

i,j=0 ' 

coll Y EmEjA = col(^2 S ^ E 



i,j=0 



N-l 



i=0 



(17) 



It follows that 



S(E 



kf: 



N-l 



k,e=o 



N-l 



N-l 



N-l 



On the other hand, 



JV-l 



J2 E^®«S(E 

k,e=o 

N-l 

(I ® Stf) (E H <g> E fc4 ) . 

i,fc/=0 



AT-1 



/, ' ) — E/,, eg) I S,/ E/,, 

fc/=0 V i=0 



(18) 



T = ^ coi(S«)coi T (E tf ) = (l®S, £ )coZ(I)coZ T (E rf ) 



i,£=0 



i,£=0 
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JV-1 /N-l \ N-l 

= i 1 ® s «) efc ® efc )W ® e * T ) = S (I ® S tf ) (E M ® E w ) (19) 

i,i=0 \k = J i,e,k=0 

QED 

In the form coZ(S/^)coZ T (E^) the derived supermatrix T appears to have first 
been studied in connection with superoperators by Jordan and Sudarshan ]TJ|, whereas the 



form [5(Ey] ^-^ was first used to give an intrinsic characterization of completely positive 
superoperators by Choi For this reason we shall henceforth denote it by Choi(S) = 
^fj=o(Ejj <8> £g> Ejj). The next Lemma will enable us to show that in the cases of 

interest here, it is a Hermitian matrix. 

Lemma 5 A superoperator S commutes with the operation of taking its adjoint, i.e. S(Z^) = 
(S(Z)) for all operators Z in its domain, if and only if it maps self-adjoint operators to 
self-adjoint operators, and if and only if for any matrix representation S of S, 

S = KSK , (20) 

where the overbar denotes the complex conjugate and 7C is defined as in Eq. jj^). 

Proof. Clearly if S commutes with the adjoint, it maps self-adjoint operators to the same. 
Now suppose that S is a matrix representation of S, and let X £ ^^xtv sa ^j s fy x = X^; 
then 

coZ(X) = Kcol(X T ) = Kcol(K) (21) 

and hence if S preserves Hermiticity, 

ScoZ(X) = coZ(Y) = KcoliY) 

(22) 

= KScol(X) = 7CS7C coZ(X) . 

Letting X range over any Hermitian basis of <C NxN now proves Eq. (pOl). And finally, if 
Z G C NxN is any (not necessarily Hermitian) matrix and «S satisfies Eq. (|20|) , we have 

col(S(Z^)) = 5coZ(Z t ) = SKcol(Z) 

K ' _ (23) 

= KScol{Z) = coZ((<S(Z)) f ) 

which, since it holds for any representation «S and matrix Z, proves <S(Zt) = {S{Z)) ] . QED 

Corollary 6 If a superoperator S commutes with the adjoint operation on its domain, then 
any Choi matrix for it is Hermitian. 



9 



Proof. Clearly a Choi matrix [<S(Ey-)] _ Q is Hermitian if and only if <S(Ey) = ((S(Ejj))^ 
for all < i, j < N, and if «S is the corresponding matrix representation of S, Lemma || 
implies 

col((S(E ji )) 1; ) = KcotiSCEji)) = KScoliEji) 

V 1 _ (24) 

= (KSJC) col(Eij) = ScoliEij) = coi(5(E y )) . 

QED 

Theorem 7 (Choi [0) Let S be a superoperator which commutes with the adjoint operation 
on its domain. Then S is completely positive if and only if the Choi matrix associated with 
any matrix representation of S is positive semidefinite. 

Proof. Let «S be a matrix representation of S and T = Choi(S) be its Choi matrix. This 
is Hermitian by Corollary and accordingly, we let 

N 2 -l 

T = WSWt = ZnU n ul (25) 

n=0 

be its eigenvector decomposition, where ti is unitary and the eigenvalues £ n are real. Then 
if £n ^ for < n < N 2 — 1, we let T n be the sequence of N x N matrices such that 
col(T n ) = y^u n . It now follows from Proposition ||]that 

(N 2 -l \ / N 2 -l \ 

T n ®T n )coI(p) = collYl T nP T i • (26) 
n=0 / \ n=0 / 

The right-hand side provides a Kraus operator sum representation for S, which by the 
previously mentioned work of Kraus [|IIJ] proves that S is completely positive, as claimed. 
Conversely, if S is completely positive, it may be expressed in Kraus operator sum form 

as 

M 



S{p) = Y, S mPSl, (27) 

m=0 

and it follows from Eq. (0) that any matrix representation «S thereof satisfies 

M 

S = Sm ® S m (28) 



m=0 
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for suitable S m G C NxN . By Proposition |3], therefore, corresponding Choi matrix T is a 
sum of dyads, that is 

M 

T = col(S m )col\S m ) , (29) 

m=0 

which is necessarily positive semidefmite. QED 

Corollary 8 Any Kraus operator sum S(p) = Xlm=o P can ^ e written in canonical 
operator sum form as 

N 2 -l 

S(P) = T nP T n > (30) 

n=0 

with (T n ,T n> ) = tv(T^T n ,) = for all < n ^ n' < N 2 and \\T n \\ 2 = (T n ,T n ) = for all 
n > M. Subject to this condition, the canonical form is unique up to the overall phase of 
the operators T n unless the Hilbert- Schmidt norms satisfy \\T n \\ = \\T n /\\ for some n' ^ n, 
in which case it is only unique up to unitary linear combinations of the operators in such 
degenerate subspaces. 

Proof. Implicit in the proof of Theorem 0. QED 



IV. A CANONICAL LINDBLAD REPRESENTATION 

We now turn our attention specifically to quantum dynamical semigroups, which (as 
mentioned in the Introduction) may be assumed to be given in the form of a superoperator 
exponential S = Exp(—J r t). The time-independent generator will usually be of the form 
T = Q + iH for superoperators Q and H, where i 2 = —1, H(p) = [p, if] for the 
Hamiltonian H of the system in question, and Q is known as the relaxation superoperator 
]Tj| . Although Q may often be self-adjoint, this is not necessarily the case. 

An important property of physically meaningful operations on density operators p, which 
we have neglected up to now, is that they preserve the trace tr(p) = 1. Given an opera- 
tor sum representation S(p) = J2m n=o Tmn P^n > ^is is easily seen to be equivalent to 
Y^m n=o Tmn Tn T m = I, the identity. We seek an equivalent condition in terms of a given 
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matrix representation S. To this end we expand S versus the basis of elementary matrices 

as 

N-l N-l 

S = Y 4( E «® E *0 = Y ^(e^e^ie^e^ (31) 

i,j,k,e=o i,j,k,e=o 

where 

4 = tr((E,,®E fcl ) t 5) = (e,®e fc ) t 5(e,®e i ) , (32) 
so that the corresponding operator sum representation becomes 

Scol(p) = coll s m e *P e #) • ( 33 ) 

For future reference, we note further that the associated Choi matrix is given by 

N-l N-l 



Choi(S) = J2 Y s% (E mn E ej (g)E kt E mn ) 

i,j,k,£=0 m,n=0 
N-l N-l 

= Y s ke (~ E ij® E ke) = Y 4 ( e » ® e fc) ( e i ® e ^) f 

i,j,k,e=o i,j,k,e=o 



(34) 



This shows that while the representative supermatrix S in this basis is formed by identically 
ordering the upper and lower index pairs of s^ e and using the resulting list as the row and 
column indices, the Choi matrix is obtained by ordering the right and left index pairs and 
using the result as the row and column indices, respectively. 



Lemma 9 A superoperator S with representative matrix S = [ s ke\ke-ij-o versus a Hilbert 
space basis {ej}^ 1 preserves the trace of its operands if and only if 

N-l 

col\l)S{p i ®e j ) = Y s ^k = Sij forO<t,j<N, (35) 

fc=0 

where 5' 1 ^ is a Kronecker delta. 

Proof. The usual trace-preservation condition can be written as 

N-l N-l 

1 = Y S M E i £ Eki = Y s kk E i* 

i,j,k,i=0 i,j,k=0 

N-l /N-l \ 

= Y [Y (efc®e fc ) f «S(e, ® O (36) 

i,j=0 \ k=0 J 
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N-l 

= (co/ t (I)5(e i ®e J ))Ey , 

i,j=0 

which is equivalent to the stated conditions. QED 

The Lemma can be stated more succinctly by saying that co/(I) is a left-eigenvector of 
S with eigenvalue 1. We note that for another important class of superoperators, namely 
the identity preserving or unital superoperators, the operator sum representations satisfy 
I = S(I) = EmnVmti^], may also be characterized in terms of their supermatrix 
representations by X^o 1 s te = ( e & ® e e)^ S col (I) = 5ke, i.e. co/(I) is a right-eigenvector of 
«S with eigenvalue 1. If S = S(t) is a unital QDS, it is easily seen that the corresponding 
Lindblad operators must be normal (or commute with the adjoints). 

Returning now to the problem of deriving a Lindblad representation for a QDS S(t) from 
a matrix exponential representation S(t) = Exp^—Ft) thereof, the obvious way to proceed, 
given the results of the previous section, is to simply differentiate it: 

JV-l 

d t Exp{-Ft)\ t=Q = -T = Yl yg(E«®E w ) 

i,j,k,e=o 

N-l N-l W 

= d t 4( E «® E w)| = 4( E ^i® E ^) 

i,j,k,i=0 i,j,k,£=0 

(note that the generator is actually time-independent). Differentiation of our trace- 
preservation condition similarly yields J2k=o $kk = ~J2k=ofkk = f° r a ^ < i, j < iV, 
and hence 

N-l N-l 

i,j,k,£=0 i,j,k,£=0 
N-l 



= Pit) ~ \ £ 4(%P + P%) (38) 

i,j,k=0 

N-l 

= E **ke ( E fc« P(0 E ^ — |(Ej£E fc jp + pE je E ki fj . 

i,j,k,e=o 

Thus we could in principle obtain a canonical Lindblad representation for p{t) simply 
by diagonalizing the (time-independent) Choi matrix of the generator [s ke ] ^ k .j\_ = 
=o fn u l™ u in] ik . je=0 , an d letting the Lindblad operators be defined by the matrices 

N-l 

i,k=0 
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- providing that the eigenvalues y?™ > for all < m, n < N. But then our trace- 
preservation condition for p implies 

JV-l N-1 N-1 

E ( L ™) tL « = E < E u Z<E je E ki 

m,n=0 m,n=0 i,j,k,£=0 , , 

JV-l N-1 / \ N-1 ^ ' 

= E E krcc ^ = E 4% = o, 

i j,t=0 m,n=0 \ / ij',fc=0 

which contradicts the fact that a nontrivial sum of positive semidefmite matrices cannot 
vanish. This result is easily shown to be independent of the choice of matrix basis. 

It follows that there must be some redundance in our choice of coefficients in any nondi- 
agonal Lindblad-type equation of the form given in Eq. fl38l) . Moreover, such an equation, 
by its very form, is assured of preserving the trace (0 t tr(p(£)) = tr(p(£)) = 0), so that the 
trace-preservation condition satisfied by the derivatives of the coefficients in an operator 
sum representation is not needed. Our problem is to find a way to modify the matrix of co- 
efficients [s^] J_ Q , while preserving the underlying mapping p i— > p, such that the result is 
positive semidefinite and so can be diagonalized to obtain a canonical Lindbladian. Because 
any Lindblad operator of the form L = a I with a G C adds nothing to p, we shall seek to 
eliminate the corresponding degree of freedom from the coefficients. 

Lemma 10 In any quantum dynamical semigroup with exponential representation S(t) = 
Exp{—Tt), the generator's matrix J- versus a Hilbert space basis satisfies 

col\l) Choi{-J=) co/(I) < 0. (41) 

// the generator is of the form T = Q + iH where Ji is a commutation superoperator and 
\G,C) = tr(C^C) = for any commutation superoperator C , then the corresponding matrix 
projection satisfies 

V l Choi(Q)V l = V Y Choi(T)V Y (V 1 = I ® I - coZ(I) coI + (I)/JV) . (42) 

Proof. To prove Eq. (f|l|), we first observe that 

JV-l / JV-l \ 

(I) Choi(-F) coZ(I) = - E (e m ®e m ) t E fki ( E y ® E kt) ) ( e n ® e n ) 

m,n=0 \i,j,k,£=0 / , , 

JV-l / JV-l \ JV-l ^ ' 

= ~£ E (/2((e| n %e n )®(et l E w e n ))) = - £ C = " tr ^ ■ 

m,n=0 \ i,j,k,£=0 / m,n=0 
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Since Choi(—T) is Hermitian, this quantity is real, and since S{t) is bounded, the eigen- 
values of — J- must all have negative real parts, so that — tr(J 7 ) < 0. 

To prove Eq. (K), we first note that it is sufficient to prove this for the commutation su- 
peroperator of an arbitrary elementary matrix E^, and transform its generating supermatrix 
to the corresponding Choi matrix: 

N-l 

Choifiij ® I - I ® Eji) = £ (( E ^ E *i) ® E ^ _ E ^ ® ( E ji E ^)) 



fc,£=0 
2V-1 



(44) 



2j (E fe j g) E fc i — E ifc <E> Ej fc ) 



fc=0 

Plugging the second term into the projection now yields 

(N-l \ AT-1 r AT-1 

£ E lk ® E jfc P 1 = E 4fc ® E ifc - coi(I) J] (e fc ® e fc ) f 
fc=o / fc=o ~ fc=o (45) 

- (e, <g> e,) co/^I) + 6 jt col(l) col\l) . 

Since the first and third terms as well as the second and fourth terms on the right-hand side 
differ only in sign, this projection vanishes identically. A similar calculation shows that the 
projection of the first term on the right-hand side of Eq. ([44]) likewise vanishes, establishing 
the Lemma. QED 

Henceforth, we take Q = T — ihi where TC is the commutator part of J 7 , and let g k \ be 
the corresponding array of coefficients. A final technical Lemma will be needed to prove the 
first real result in this section. 

Lemma 11 // (g,C) = ti&C) = for every commutation superoperator C as above, then 
the coefficients g\\ of any supermatrix representation Q satisfy 

N-l N-l 

£ g& = E 9$ (46) 

k=0 1=0 

for all < m,n < N. 

Proof. The proof is by direct computation: 

(N-l 
(E nm <g> I - I <g> E mn ) 9k% ( E ^' ® E feO 

i,j,k,i=0 
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tJV-1 
all ((Enm E Ej ) ® E fci - E# ® (E mn Efci)) 

(N-l \ / JV-1 

E ^ ( E ^- ® Efc >) = tr E ^ ( E ^- ® E ™) 

i,j,k=0 / \ i,j,i=0 

jV-1 iV-1 

E tr ( E ^) tr ( E ^) = E tr ( E «) tr(E mi ) 

i,j,k=o i,j,e=o 

N-l N-l 

E g km $nj 5ki = E 9ne $tj $mi 

i,j,k=0 i,j/=0 



(47) 



QED 



This Lemma may be paraphrased by saying that the "partial trace" (or contraction) 
of Q with respect to either its left or right Kronecker factors are the transposes of one 
another. 



Proposition 12 Let S(t) = Exp(-J-'t) be a quantum dynamical semigroup with T — Q + 

N-l 



TV- 1 



iH as above. Then if their supermatrices versus a Hilbert space basis are J- = \fu\kt-ij- 
G = [9k£\ki-ij = o an< ^ ^ = I®H — H(g> I, we have 



N-l N-l 

Pit) = E fSEupWEjt = i[p(t), H] - E 9%EkiP{t)E je 

i,j,k,e=o i,j,k,e=o 

N-l 



[p(t), H] - - E 9u (2Efeip(t)E^ - E je E kl p(t) - p(t)E#E w ) (48) 

JV-1 / iV-1 1 JV-1 

[p(t), H] - E E ^E fei p(t)E„ - -E^S (e^pW + p(t) 



Eji 



i,j=0 \ k,£=0 k=0 

N-l 



where \9ht\ki-ij = o are ^ e coe J^ c ^ en ^ s of the supermatrix 

Q = Choi(V l Choi(Q) V 1 ) = Choi(V l Choi(F) V 1 ) . (49) 



Proof. Note that H occurs on both sides of Eq. ( f48"D , so we can just ignore it (i.e. set H = 0) 
in the proof. Since col(T)col\T) = Y^l~n=o E m« ® E mn , we find that: V 1 Choi(Q) V 1 

N-l / ^ N-l 

= 9u I Ejj <8> E^ — — (E mn Ejj (g E mn E^ + EjjE mri <g> E^E mn ) 



i,j,k,l=0 \ m,n=0 



1 iV_1 \ 

^ ^ E mn EjjEpq (gi E mn E^£Epq I 

m,n,p,q=0 / 
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JV-1 / ^ / JV-1 JV-1 x 

gH ®ij ® E fc£ - — f 5* ^ E mj <g> E m ^ + ^ ^ E m <g> E kn J 

i,j,k,i=0 V ^ m=0 n=0 ' 



(50) 



A-/ 



m,n=0 

= E U ■ | E -C - | E A + $ E .»=) % ®e 

i,j,k,£=0 ^ m=0 n=0 m,n=0 ^ 

JV-1 

= E ^%®E M . 

i,j,k,£=0 

Equation flSOD thus determines the projected coefficients i n terms of the original coeffi- 
cients, and if we compute the Lindbladian versus a Hilbert space basis using the projected 
coefficients as in the last line of Eq. (f48D, we get 

JV-1 / JV-1 JV-1 \ 

C (P) = E E 3« E« p(t) - - J] S« (E rt p(t) + p(t) B i4 ) . (51) 

i,j=0 \k,£=0 k=0 / 

The supermatrix representation of the first operator sum in this equation can be further 
simplified as follows: 

N-l N-l N-l / „ JV-1 



XV— X IV— X iV — 1 / ^ — 1 \ 

E 0& E «® E w = ^E £i ®E fei - ^ (_^^W<g)I 

„,j,M=0 i,i,fe,i=0 j,£=0 ^ m=0 ' 

jv-i / jv-i v / jv-i v y ) 

- E hvE^)i®E fei + 1 — «c)i®i- 

i,fc=0 ^ n=0 ' ^ m,n=0 ' 

Similarly, by Eq. (|50|) the supermatrix representation of the second operator sum in Eq. (|5T|) 
simplifies to 

2 E s k*$kt + I®E y ) 

ij,fc/=0 

JV-1 /jv-i 1 jv-i 1 JV-1 \ 

- 2 E E - n E ^ - n E (E^I + I®^) (53) 

i,j=0 \ k=0 m=0 n=0 / 



/ 1 \ 

+ [jp E ffZji®i 



m,n=0 

Taking into account the difference in the signs of the operator sums in Eq. (|51|) , the last 
terms on the right-hand sides of Eqs. (j52|) and (^) clearly cancel, while the first summation 
on the right-hand Eq. (|53|) vanishes by our trace preservation condition (cf. Lemma |9|). The 
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remaining terms on the right-hand side of Eq. (|53| ) can be rearranged using the symmetries 



of the summations proved in Lemma [11], as follows: 

JV-l /JV-l N-1 

2N 



^ E E & + £ C ) (E* + e,o 

j j=0 \ m=0 n=0 / 

JV-1 / .JV-l v .JV-1 , 

iV E E ^ E 0- ® I + £ & (54) 



2iV 

* J = V m=0 m=0 

+ E^nT E^n 
^ n=0 ' ^ n=0 ' 

.v-j , .JV-1 v .JV-l v 

y E £C + E4 n )i®e 



i,.7=0 \ N m=0 n=0 

It is now apparent that these terms cancel with the second and third terms in Eq. ( |52| ) after 
a change of dummy indices, leaving only its first term behind. QED 

Thus, roughly speaking, the subtraction of the (L^Lp + pL^L)/2 terms from the LpU 
terms of the Lindbladian simply ensures the trace of p still vanishes after projecting out the 
commutator and identity superoperator parts of the derivative of the corresponding Kraus 
operator sum. It remains to be shown that the Choi matrix of the operator sum is positive 
semidefinite if and only if the Choi matrix of the projection of its derivative is positive 
semi definite. For the sake of completeness, we first prove the following (well-known) result, 
using only the techniques developed above. 

Lemma 13 The composition of two completely positive superoperators Ao B is again com- 
pletely positive. 

Proof. Let W Diag(at) Vl^ and VDiag(/3) be the eigenvalue decompositions of the 
supermatrices Choi(A) and Choi{3) respectively, and consider the Choi matrix of their 
product, namely 

JV-l / JV-1 \ 

Choi(AB) = £ £ 

i,j,k,i=0 \ m,n=0 ) 
N-1 / JV-l . JV-l \ / JV-l 



l\ — l / l\ — 1 . — 1 \ . 1\ — 1 \ 

£ £ £ E «xK®e* 

i,k,£=0 \m,n=0 ^ p,q=0 ' ^ r,s=0 ' 

N-1 / JV-l .JV-l v .JV-l s \ 

E E E«0 E% P < JEy®E w (55) 

i,j,k,£=0 \ p,q,r,s=0 ^ m=0 ' ^ n=0 ' ' 
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N-1 /N-1 ,N-1 x \ /N-1 ,N-1 

E <W E E ^ (e,®e,)t Y, 

,q,r,s=0 \i,k=0 V m= n / / V „• />— n V ™=n 
iV-1 

E <& r <«) f 



v np v jr 

a lq U ns 



p,q,r,s=0 \ i,k=0 x m=0 / \ 3,1=0 n=0 

iV-1 



p,g,r,s=0 

Such a sum of positive semidefinite matrices (Hermitian dyads, in this case) is always again 
positive semidefinite, proving the Lemma. QED 

Theorem 14 The integral of a Lindbladian yields a quantum dynamical semigroup, and 
conversely, and the derivative of any quantum dynamical semigroup can be placed in canon- 
ical Lindblad form. 

Proof. Given any Kraus operator sum for a quantum dynamical semigroup S(t), we know 
that its time-derivative will be equal to the result of applying a fixed generator —T to the 
density operator p(t) at any given t > 0. Integration of a matrix representation thus yields 
S(t) = Exp(-^Ft), and for a sufficiently small 5t > this exponential may be approximated 
arbitrarily closely by 

Exp(-Ft) w X - J=bt + 0((5t) 2 ) = X - Q8t - iHSt + 0((5t) 2 ) , (56) 

where X = I <g) I and iH. denotes the commutator part of J- . Since S(t) is completely 
positive, any Choi matrix for it must be positive semidefinite, and so must any projection 
thereof, in particular, 

V x Choi{X - T bt)V x = -V I Choi(g)V 1 5t = - Choi{Q) St . (57) 

This allows — Q and hence also its sum with — iH to be placed in canonical Lindblad form, 



which by Proposition must have the same action on any p as the differential superoperator 
-JF. 

Conversely, suppose that a given superoperator T — Q + iH can be placed in canonical 
Lindblad form, 

M 

-HP) = C(p) = -»W(p) + E (L mP Ll-\LlL m p-\pLlL m ) (58) 

m=l 
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where TC(p) = [H, p] for the commutator part of T . In terms of a matrix representation H, 
L m of these operators, this is equivalent to 

M 

_jp = -g- lH = ^(L^^-il^L^-iLX^l) (59) 

m=l 

+ «(H®I-I®H) . 

Then over a sufficiently small time interval St, the exponential (integral) can be approximated 
arbitrarily closely by the product of the exponentials 

, M v 

^ m=l ' 



Seep ^ L m ® L^j Exp(—8t itL ) + 0((5t) 2 ) 



m=l 



(60) 



= A(5t) B(5t) C(5t) + 0({5t) 2 ) . 

Since the two types of terms in the argument to the first exponential commute, it evaluates 
to a Kronecker product, namely 

M 



m=l ' 

( M \ ( M — — \ 

= Expi- \8t J2 L m L m) ® ^P(- §** E 

m=l ' m=l ' 

Thus by Lemma [| the corresponding Choi matrix is the dyad 

Choi(A(8t)) = col(Exp(-\8tY,^^) col^Expf-lSt^^r 

^ ^ m=l ' ' ^ ^ m=l 



(61) 



(62) 



which is necessarily positive semidefinite, proving that A(t) is a QSD all by itself. As for 



the second factor in Eq. |6(| we may expand it as 

M \ M 



5t^2 L m ®L m ) « I + L m ® L m + 0((8t) 2 ) 

m=l ' m=l 



(63) 



Because V 1 Choi(X <g> I) V 1 = V 1 Choi(I ® X) P 1 = for all X G C NxN , the Choi matrix 
of the summation on the right-hand side is easily seen to be Choi(Q) = "P 1 Choi(Q) 7* , 
so that 

Choi{B{5t)) w 00^(1)00^(1) + 5tChoi{g) + 0((5t) 2 ) . (64) 
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The Choi matrix Choi(Q) = ^ ? ^ =1 co/(L m ) co/^(L m ) is of course positive semidefinite, 
and (since Q = J^ =1 7n(U n ® U n ) where j n , co/(U n ) are the eigenvalues and eigenvec- 
tors of Choi(Q)) so are the Choi matrices of all higher terms in the Taylor expansion of 
Choi(B(5t)), thus showing that £>(£) is also completely positive for all t > 0. Finally, the 
last factor of Eq. ( |6Qj) , 

C(5t) = Exp(-dttn) , (65) 

is unitary and hence likewise corresponds to a completely positive superoperator for all time. 

It now follows from Lemma [13| that for 5t <C the product of all three factors 

A.(5t), B(8t), C(5t) in Eq. (|60D will be completely positive, and hence for any given t > 
the telescoping product 

Exp(-tF) « ( Ajt/n) B(t/n) C(t/n) ■ ■ ■ Ajt/n) B{t/n) C(t/n\ ) + 0{{t/nf) (66) 

n times 

will also be completely positive for all n > t/5t. The Theorem now follows by noting that 
the set of completely positive superoperators is closed, and taking the limit as n — >• oo. QED 



V. APPLICATION TO QUANTUM PROCESS TOMOGRAPHY 

The applicability of the foregoing results to QPT derives from the following theorem, 
whose origins can be traced back to work by Eckart, Young and Householder JL7|, [18) , 
and has since given rise to a field of statistical data analysis widely known as "principal 
component analysis" flI5| . The present author has proven it several times in the course of 
his career |2U], ^T], [EJ, and regards the following proof as the simplest. 

Theorem 15 Let M e C AfxW be a Hermitian matrix with eigenvalue decomposition 

N-l 

M = U f AU = J2 X i u i u U ( 67 ) 
e = o 

where the eigenvalues have been sorted in nonincreasing order Xi > A^ + i for < i < N — 2. 
Also let P denote the convex cone of positive semidefinite matrices in <C NxN and 7-p(M) be 
the orthogonal projection of M onto P with respect to the Hilbert- Schmidt (or Frobenius) 



matrix norm \\ ■ \\, which satisfies the "least- squares" criterion jf3|/ 



M - VJM) II 2 = min II M - M' f . (68) 
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Then we have 

N*-l 

V F (M) = M* = U f A* U = \*u £ ul , (69) 

e = o 

where A* is the diagonal matrix of eigenvalues A with all of its N — N* negative eigenvalues 
set to zero. 

Proof. Any positive semidefinite N x N matrix can be written as XX^, where X G C NxN ' 
and N' is its rank. It follows that the minimum in Eq. (p5|) can also be written as 



min C(X) = min II XX 1 " - M f . (70) 

XgC iVxJV' Tt£C NxN ' 

It is easily seen that the gradient matrix of ((X) is 

^ = ^tr((XXt - M) 2 ) = 2 (XXt - M) X . (71) 

On setting this to the zero matrix, we obtain the nonlinear matrix equation 

MX = X(X f X) , (72) 

wherein X^X is an N' x N' Hermitian matrix which, for N' = 3 and X G M Arx3 , is essentially 
the inertial tensor (plus a multiple of the identity) of a system of unit mass points with 
coordinates e^X (0 < i < N). Since the Hilbert- Schmidt norm is unitarily invariant, we 
may assume that these "coordinates" have been chosen so that X^X = Diag(£i, . . . , £at/) is 
diagonal, in which case Eq. ([72]) becomes 



Mx, = (j = 0,... ,N-1), (73) 

where Xj = Xe.,- are the columns of X. It follows that the Xj are proportional to the 
eigenvectors Uj associated with certain nonnegative eigenvalues Xj = £j of M where, since 
||xj|| 2 = Q , the constant of proportionality is v/Aj ■ On expanding the trace in the function 
C(X), we now obtain 

C(X) = tr(M 2 - 2XX f M + (XX f ) 2 ) 

= tr(M 2 ) - tr(2X t MX - (X f X) 2 ) . 

By Eq. (|?2|), however, the matrix X G C NxN ' that minimizes C(X) satisfies 

(X) f MX = ((X) f X') 2 = Diag(A 2 ,... , A 2 ^ ) , (75) 
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(74) 



so that 

JV'-l 

C(X) = tr(M 2 ) - X l ( 76 ) 

3=0 

From this we see that, for any integer N" with < A^" < N' and Xj > for < j < N", 
the minimizing X" G C NxN " is obtained by setting X" = [s/X^Uj] =Q . It follows that the 
minimizing positive semidefinite matrix X* (X*)t is obtained by setting N" to the number 
N* of positive eigenvalues of M. QED 

This theorem can be used to "filter" statistical estimates of either superoperators or their 
generators so as to obtain a completely positive estimate. In the case of an estimate S' of a 
matrix representing an unknown superoperator S, one simply sets any negative eigenvalues 
of the associated Choi matrix T' = Choi(S') to zero, reconstructs the improved estimate 
T* from these modified eigenvalues and the original eigenvectors as in the theorem, and 
converts the result back into a new estimate «S* = Choi(T*) of the superoperator via the 
same involutory mapping Choi. The theorem assures us that this procedure makes the 
smallest possible change in T', with respect to the Hilbet-Schmidt norm, so as to render it 
positive semidefinite and so ensure that «S* represents a completely positive superoperator. 
Because the mapping Choi simply permutes the elements of its argument, we can be sure 
that this procedure also minimizes the change \\S' — S*\\ in S'. We now show that <S* is 
assured of being an improved estimate of the corresponding matrix of the true superoperator 
S, again in the least-squares sense. 

Corollary 16 For S, S' and S* defined as above, we have 

|| S* - S || < || S' - S || . (77) 

Proof. Since S* is the orthogonal projection of S' onto the convex cone Choi(¥) of ma- 
trices representing completely positive superoperators, S 1 — S* is orthogonal to a supporting 
hyperplane at S*, while by its definition S G Choi(F) must be on the opposite side of this 
hyperplane from «S'. This in turn implies that the angle 6 between S and S' at S* satisfies 
9 > 7r/2, and hence by the law of cosines 

> cos(#) = idlS^-S'f + ||S*-S|| 2 - \\S'-S\\ 2 ) , (78) 

i.e. \\S' - S\\ 2 > \\S* - S'\\ 2 + ||S* - 5|| 2 > \\S* - S\\ 2 . QED 
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The procedure in the case of a QDS generator T = Q + iH, is a bit more involved, since 
one needs to compute the projection of the Choi matrix £' = "P 1 Choi(J-') "P 1 of the esti- 
mate J-' before diagonalizing it. This of course will remove the Hamiltonian superoperator 
component, which must then be obtained by some other means. In addition, one cannot 
reconstruct a matrix Q* for the decoherent component Q of T from the matrix £* obtained 
by setting any negative eigenvalues e m of £' to zero simply by applying the Choi mapping, 
since the other terms needed to preserve the trace will also have been lost in the projection 
(if indeed the estimate J 7 ' itself were trace-preserving). Instead, one has to construct all 
the Lindblad operators L m such that coZ(L m ) = y/e^,v m , where e m > 0, v m are eigenvalue, 
eigenvector pairs of £', and compute Q* as indicated in Eq. (|6]). As a result, there is no 
guarantee that Q* will be closer to its true value Q versus the Hilbert-Schmidt norm, al- 
though we expect that this will usually be the case. Further discussion regarding how one 
might go about solving these problems must take the exact experimental situation at hand 
into account, and as such is outside the scope of this paper. 

In the remainder of this section we will illustrate how the above results may be applied to 
a simple example, namely the Block equations for a single spin 1/2 qubit in a frame rotating 



at its Larmour frequency in an applied magnetic field [[HJ. As is well-known [23], these can 
be expressed in canonical Lindblad form as 

p = £(p) = (2EoipE 10 -Eoo/o-pEoo) + 

(2E 10 pE i-E 11 p-pE 11 ) + (79) 



l-A 
4Ti 

\ 2T 2 4Ti 



((Eqo-Eu)p(Eoo-Eu) - P ) , 



where T\ and T2 are the characteristic relaxation and decoherence times and A = po — p± 
is the excess probability in the ground state E o at equilibrium. The supermatrix of the 
generator versus a Hilbert space basis in the ordering E o , E 10 , E i , En induced by the 
"coZ" operator is 



l+A 

4Ti 

l-A 

4Ti 



(2 E01 ® E01 - I ® E 00 - E 00 <8> I) + 
(2 E10 ® Em - I ® E n - E n <g> I) + 



{m--m){ ( E °° " El1 ) ® ( E °° - En) - 1 ® : ) (so) 
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J-A Q Q 1+A 

-A 





1-A 



~% 





2Ti 




1+A 



2Ti " " 2Ti J 

The time-dependent exponential of this matrix may be shown to be 



Exp( - Ct) 



(l+e-'/ T i)+A(l-e- t / T i) (l- e -'/ T i ) +A (l- e -*/ T i ) 

2e"*/ T 2 

2e~*/ T 2 

(l-e-*/ T i)-A(l-e- t / T i) (l+ e -'/ T i ) - A (l- e -*/ii ) 



(81) 



which in turn corresponds to the Choi matrix 

'(l+e- i / T i)+A(l-e- t / T l) 



M(t) 












2e-'/ T 2 
(l-e-*/ T i)(l-A) 

(l-e-*/ T i)(l+A) 

(l+e- t / T i)-A(l-e- t /ri)_ 



(82) 



This in turn is readily shown to be positive semidefinite for all t > if 27\ > T2 . Its 
derivative at t — 0, however, is 



M(0) 



Choi(C) 



2Ti u 

X - A 





T 2 





2Ti 

^ 

u 2Ti 





_ J_ 




1+A 
2Ti . 



(83) 



and the outermost 2x2 block of this matrix is positive semidefinite only if 27\ < T 2 . 
Applying the projection V 1 = X — col(T) coV(T)/2 converts it to 



-V 1 Choi(C) V 1 



_2 i_ 

2T 2 4Ti 






1 

,4Ti 



1 
2T 2 





1-A 
2Ti 










1+A 
2Ti 





1_ 

4Ti 2T 2 




1_ 

2T 2 4Ti 



(84) 



which is now positive semidefinite with eigenvalue, eigenvector pairs: 
0. 



1 1 

2T 2 4Ti ' 



1+A 
2Ti ' 






\ fl-A 





1 










J ' I 2Ti ' 


1 









(85) 



The eigenvectors are easily seen to be obtained by applying the "coZ" operator to the 
matrices I = E o + En , E o — En , Ei and E i , returning us to the canonical Lindblad 
form in Eq. (|79|). 
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We will now use this example to illustrate how the matrix formulae obtained in this paper 
can be applied to QPT, by numerically simulating the "data" needed for QPT from the 
above solution to the Bloch equations. These data correspond to an experimental scenario 
in which a set of precisely known input states {p'k}k=i were allowed to evolve under the 
propagator in Eq. ( |5T]) for varying periods of time, and the results {Pfc Ut } fc=1 determined 
by state tomography 0, ||, |25], |2B], [27). Assuming that the input states span the space of 
single-qubit Hermitian operators, this allows us to determine the propagators at each time 
point according to 



Exp( -Ct) [col{p"), . . . , col(p-)] = [col( P r), • • • , col{pT) A 

(86) 

^Exp(-Ct) = [col( P r),... ,col(pT)][col(p[ n ),... ,col(p m K \ 



Although this relation is exact when the output states are known precisely, in actual prac- 
tice experimental errors would result in only an approximate estimate S'(t) of the actual 
propagator S(t) = Exp(-Ct). If one obtains such estimates at an arithmetic sequence of 
time points = to, ti, . . . , tj = Jt\, however, one may solve a linear least-squares problem 
to obtain an improved estimate of the propagator «Si = S(ti) at the first nonzero time point 
||, namely 

j-i 

min T (x(T)) where X (T)= ^ \\T - S' j+1 || 2 . (87) 

3 = 

One may of course set «Sq = S(t ) = I ® I, the 4x4 identity, and one should also filter 
the remaining estimates by symmetrizing their Choi matrices (i.e. by adding them to their 
adjoints and dividing by two), setting any negative eigenvalues ip = and transforming back 
to a new estimate (as described previously). The minimizing solution to this least-squares 
problem is easily shown to be 

s'[ = (]T ^(^(E 3(3)0*' (88) 
j=i j=i 

where in most cases the Moore-Penrose inverse (f) may be replaced by the usual matrix 



inverse 



Finally, S" may be converted into an estimate of the generator via the matrix "pseudo- 
logarithm", Plog. This is computed by diagonalizing S'[ = W$W _1 , setting any eigen- 
values <fti < or (pi > 1 to zero while taking the usual logarithm of the rest, then performing 
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the inverse similarity transformation and dividing by ti, i.e. 

3 

t x C" = Plog(S'i) = WPlog^W- 1 = ^ plog(&) (WeJ (W^e,) 1 , (89) 



i = 



where 

if < fa < 1; 

plog(&) = { (90) 

otherwise. 

The eigenvalues will be real since no Hamiltonian was assumed in the simulations, and 




arguments similar to those involved in Theorem |15| can be used to show that the pseudo- 
logarithm will then yield a generator C" that minimizes \\S'[ — Exp(—C"ti)\\. Lastly, 
the estimate C" is filtered by projecting its symmetrized Choi matrix by 7-* 1 , setting any 
eigenvalues e = 0, and reconstructing to obtain the optimum estimate £*, as described 
above. 

The specific values of the parameters used for the simulations were T\ = 0.5, T 2 = 0.1 
and A = 0.1; the relaxation times T\ and T 2 are typical of liquid-state NMR samples, while 
the polarization A was deliberately made larger to render it visible despite the noise. In 
accord with Eq. flSCP, these gave rise to the generator 



-0.9 o 1.1 
o -10.0 o o 
o -10.0 o 
0.9 -1.1 



(91) 



The input states were taken to be E o, En, (e + ei)(e + ei)/2, and (e — iei)(e + iex)/2, 
while the times used were set to tj = j/4 (j = 0, . . . , 4). Finally, the noise levels evaluated 
were Qi = 0.01, Q 2 = 0.05 and f2 3 = 0.25, where the noise was simply added to the output 
states {pfc Ut } with a Gaussian distribution, zero mean, and variances <rj Vl\ proportional to 
the mean-square size cr| of the elements of S(tj) (j = 0, . . . , J; k = 1, 2, 3). The results below 
were averaged over 100 independent estimations of the propagators at each time point, using 
different random noise for each estimation and time point, followed by filtering and fitting 
to obtain estimates of the generator, all at each of the three specified noise levels. 

Table 1 shows the average changes made to the propagator estimates upon symmetrizing 
and filtering the eigenvalues of their Choi matrices, as measured by the Hilbert-Schmidt norm 
of the difference divided by that of the true propagator, together with the average number 
< 0}| of eigenvalues set to zero in the process. It may be seen that the changes in the 
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h 




h 


u 


IW<0}| 




0.0108 


0.0121 


0.0116 


0.0127 


0.000 




0.0581 


0.0601 


0.0644 


0.0605 


0.000 


^3 


0.3062 


0.3038 


0.3074 


0.3098 


0.290 



TABLE I: Average over 100 runs of Hilbert-Schmidt norms of the changes in the propagators on 
symmetrizing and filtering the eigenvalues {V>} of their Choi matrices, divided by the norm of the 
actual propagator ||«S|| (see text); the last column shows the average number \{ip < 0}| of negative 
eigenvalues of that were set to zero. 

estimated propagators upon filtering became significant as the noise level increased, but were 
generally little more than the added noise. Negative eigenvalues were frequently encountered 
only at the highest noise level Q 3 = 0.25, however, so in fact most of these changes were due 
to the symmetrization needed to make the estimated Choi matrices Hermitian. 

Table 2 shows the average changes made to the various generator estimates computed 
(this time normalized by the norm of the true generator), together with the average numbers 
of eigenvalues set to zero in computing the pseudo-logarithm (|{0 < 0}|) and in filtering 
(\{e < 0}|). Again, few eigenvalues with incorrect signs were encountered either in comput- 
ing the pseudo-logarithm, or in symmetrizing and filtering the resulting generators. This 
means that, once again, most of the improvement was obtained via the projection "P 1 and 
subsequent reconstruction, forcing the estimated generators C* to preserve the trace (which 
the unfiltered estimates C" did not). Finally, it should be noted that the filtered generators 
C* usually came out closer to the actual solution than the unfiltered, although this was not 
invariably so. Together, these numerical results strongly support our claim that the formulae 
derived in this paper provide a powerful set of tools with which to tackle quantum process 
tomography on systems that may be aptly modeled as a quantum dynamical semigroup. 



VI. CONCLUSIONS 

In this paper we have presented formulae by which the supergenerators and superpropa- 
gators of quantum dynamical semigroups may be manipulated, placed in canonical Lindblad 
and Kraus form, and all these forms interconverted. These formulae constitute a set of tools 



28 





\\c"-c*\\ 


\\c"-c\\ 




|{0<O}| 


\{e < 0}| 




0.0077 


0.0305 


0.0300 


0.000 


0.000 




0.0634 


0.1720 


0.1676 


0.010 


0.420 


^3 


0.2971 


0.6355 


0.5553 


0.580 


0.840 



TABLE II: Average Hilbert-Schmidt distances (columns 1-3) among the estimates of the gener- 
ators divided by the norm of the actual generator and (columns 4-5) average numbers of 
eigenvalues set to zero in obtaining these estimates (see text). 



that should be particularly valuable in developing robust procedures for quantum process 
tomography and quantum channel identification IJ, using diverse forms of experimental 
data. We have illustrated one such application using data simulated from the well-known 
Bloch relaxation equations on a single spin 1/2 qubit ||16|| , which assumed that full state 
tomography versus a basis of input states could be performed. This example demonstrated 
the anticipated robustness of the procedures employed, which was the result of combining 
the eigenvalue characterizations of completely positive supergenerators and superpropaga- 
tors derived in this paper with powerful matrix approximation methods derived from the 
field of principal component analysis [[HJ . 

It should be clearly understood, nevertheless, that the procedures given here were in- 
tended primarily to provide a concrete example of how the mathematical results given in 
this paper can be applied to quantum process tomography, and not as a prescriptive recipe 
that is in all cases optimal — or even applicable. For example, the system of interest will 
often evolve coherently as it relaxes towards equilibrium, and at a rate far larger than the 
relaxation processes themselves. In this case the relaxation generator itself will be averaged, 
significantly complicating its physical interpretation, and the superpropagators determined 
from full state tomography versus an input basis set will usually have complex eigenval- 
ues. Even assuming its matrix can be fully diagonalized, the well-known ambiguity of the 
matrix logarithm with respect to the addition on arbitrary multiples of 2m onto its eigenval- 
ues will render our "pseudo-logarithm" technique inapplicable. Particularly in such cases, 
better results can be expected from nonlinear fits of the supergenerator to the superpropaga- 
tors [JTT], , but the question of whether these problems are best solved by computational 
means, experimental means, or some combination thereof, will clearly depend upon the 
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circumstances. 

There are further many other ways to represent a quantum state besides a density matrix, 



for example by a Wigner distribution [ 25f , or it may even be desirable to forgo state tomog- 
raphy altogether and to base quantum process tomography on a sequence of time-dependent 
observations which, although individually insufficient to fully determine the superoperator 
or even the system's quantum state, nevertheless do so in aggregate. Alternatively, one 
might utilize a form of indirect measurement via qubits outside of, but interacting with, the 
system of interest Jl], [16], gjj. We anticipate that many creative applications and extensions 
of the techniques introduced in this paper will be developed in the years ahead, as quantum 
information processing technologies progress towards experimental reality. 
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